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[57] ABSTRACT 

Provided is an apparatus for analysing digital radio 
transmissions which can provide an output indicating 
the condition of a digital transmission link. The appara- 
tus operates by sampling received radio signals to pro- 
duce for each sampling instant a signal or signals repre- 
sentative of the modulation state of the transmission at 
the sampling instant The apparatus includes a processor 
which can receive and measure a given number of such 
signals. The processing means is arranged to process the 
signals according to one or more stored routines and 
generate one or more parameters which indicate the 
condition of the transmission. The apparatus includes 
means for representing the samples by digitally encoded 
numbers and the processing means arranges those sam- 
ples into one of a plurality of groups according to the 
value of the number. The number of groups correspond 
to the number of modulation states. This technique 
allows the apparatus to carry out accurate analysis of 
the transmission. 

4 Claims, 10 Drawing Sheets 
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The invention will be described now by way of exam- 

ANALYSIS OF DIGITAL RADIO TRANSMISSIONS pie only with particular reference to the accompanying 

drawings. In the drawings: 

This invention relates to the analysis and testing of FIG. 1 is a schematic illustration of a digital radio 

digital radio transmissions. In particular it relates to 5 transmission system; 

apparatus which can be used to carry out such analysis. FIG. 2 is a view showing the front panel of a device 

In a digital radio transmission system digitally en- m accordance with the present invention; 

coded data is transmitted from a transmitter in the form F 10 - 3 ■» a Wock circuit diagram of a device in accor- 

of modulated carrier waves. Various modulation dance with the present invention; 

schemes are known and these allow data to be transmit- 10 4A t0 40 ^ examples of constellation dis- 

ted at varying rates. For example in 16 QAM two quad- p1 ^?' ^ . » .„ . 

rature carrier waves are transmitted and each carrier mGS ' 5 t0 10 lllustrate the operation of the present 

wave can be modulated with one of four levels. This device. . _ J . , 

arrangement provides 16 modulation states and enables t€ , ™ e P resent descnpton is concerned with apparatus 

each group of four bits in the digital data to be repre- 15 testm S or ^ ysm % t™~ A 

sented by a particular one of the 16< states. P * *ff*ilu? ^wf™*? 

a * a\_ • • . A , ^ • - , . . FIG. 1. At a transmission station 10 data, which may be 

£V I? g T^r!, *™^.««^ PCM data, on line 11 is applied to a modulating device 
are demodulated and applied to decision circuitry u The modulati devi £ 12 shown m mG \ ^ de _ 

which recognises the particular moduktion states and 2Q ^ d to provide 16 QAM modulation and has two 

outputs the appropriate digits so that the digital data is paths l ^ Q Qne of which . ^ a sing waye md the 

reconstituted. other a cosme wave . Each path includes a four level 

It is important that such radio links can be analysed mo dulator 15, 16 and this arrangement provides 16 
and tested since if the transmission degrades errors can possible modulation combinations for representing the 
occur. One way of observing errors is to use the two 2 5 digitally encoded data. The modulated signals are fil- 
demodulated signals to form what are known as eye tered m fii ters 13, 14 and transmitted from a transmitter 
diagrams, for example on a cathode ray tube display. ig which may include a travelling wave tube. 

Another way in which the transmission can be ob- At a receiver 20 the transmitted signals are detected 
served is to sample the received demodulated signals at and filtered and then demodulated at demodulators 21, 
instants corresponding to decision instants and to repre- 30 22 by signals corresponding to the sine and cosine carri- 
sent the results on a display in an X-Y co-ordinate for- ers. These signals are then filtered in filters 23, 24 and 
mat. In the case of 16 QAM this produces a 4 by 4 array applied to a decision circuit 25 to reconstitute the digital 
of clusters of points, each cluster corresponding to a data which is then output on line 26. There are known 
particular modulation state. Ideally there should be 16 ways of analysing digital radio transmissions. One of 
dots on the display but in practice these appear as small 35 these is to observe the eye diagram. Another is to con- 
clusters or groupings of dots. In a system which is func- struct a constellation which is an array of clusters of 
turning correctly the display will comprise tightly points on for example a cathode ray tube screen each of 
grouped clusters of dots in a perfectly square configura- which corresponds to one of the modulations of the 
tion sometimes known as a constellation. Any distortion transmitted radio signals. In the case of 16 QAM an 
from this arrangement are due to faults and the particu- 40 array of the type shown on the display 30 (FIG. 2) of 
lar form of the distortion can be used to identify the the present device is obtained. The display is formed by 
type of fault sampling the recovered I and Q signals at a plurality of 

The present invention is concerned with apparatus instants to form a large number of samples. These sam- 
which can produce both eye and constellation displays P !ed values m displayed on the cathode ray tube dis- 
and which has a graticule selection facility together 45 P la > 30 M shown m FIa 2 can ^ be stored m 
with a simple means for adjusting the constellation to form for further analysis. The array shown is for 

the graticule. The apparatus can also provide quantita- 16 Q^M and comprises ^6 clusters 31, one for each 
tive analysis of the constellation which has not hitherto ^dulation state or combination. Each sample corn- 
been possible pnses an I value and a Q value. The I values extend 

According"to the present invention there is provided 50 ^ the x ™* md the <* v ^ ues * e Y Each 
apparatus for analysing digital radio transmissions com- P ™ tm f a clus f ter corresponds to a sample, 
prying means for sampling received radio signals to ]^\ tyve <? f f™* ™ [ b *™ d t0 .^Vse the digital 
, r ■ * * • * , radio transmission. If the transmission is functioning 

produce for each samphng instan a si^al or signals CQrrectl then a ^ non .di storted will ^ 

representative of the modulation state of the transmis- 55 generated . Faults m th e transmission such as for exam- 

sion at the samphng mstant, and processing means le recdver lock le error or transmitterA 

which can receive and measure a given number of said receiver quadrature ^ error them selves as 

sisals, said processing means being operable to process distortions in the constellation. The apparatus to be 

said signals in accordance with one or more stored described can not only display the constellation but can 

routines to generate one or more parameters which are ^so carry out measurements on the received digital 

indicative of the condition of the transmission, charac- rad i 0 transmissions to provide a quantitative indication 

terised in that said samples are represented by digitally 0 f the magnitude of errors which occur in the transmis- 

encoded numbers and said processing means is arranged sion. 

to allocate said numbers to one of a plurality of groups, The device is illustrated as a block 34 in FIG. 1 and is 
the number of groups corresponding to the number of 65 coupled to the receiving equipment between the filters 

modulation states of the transmission. 23, 24 in the I and Q paths and the decision circuitry 25. 

Cross-reference to U.S. Pat. No. 4,825,449 should be When connected in this way the device receives the 

noted. demodulated filtered I and Q channels. The device has 
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various modes which enables the eye diagram for the I operable controls coupled to the circuit 92. For each 
channel to be displayed, the eye diagram for the Q sample there is produced at the outputs of the circuits 54 
channel to be displayed, a constellation of the type and 60 an I level signal and a Q level signal. These 
shown in FIG. 2 to be displayed, and a measure mode in signals can be fed via the circuits 56 and 61 to produce 
which measurements of the samples making up the 5 a display dot on the cathode ray tube 58 at a position 
constellation can be carried out. These , measurements corresponding to those I and Q levels. By producing 
can be output either to the display 30 or to other exter- several hundred or thousand samples it is possible to 
nal equipment Details of the measurements will be generate a display of the type shown in FIG. 2. This is 
given later. The device also has a print mode which a constellation display. 

outputs constellation data together with measurement 10 The signal levels of each sample can also be fed via 
results to an external printer. The mode can be selected the re-sampler 65 and digitised using the comparator 
by operation of the buttons shown in the region 35 of which forms part of that circuit and a digital to ana- 
the panel in FIG. 2. logue converter which forms part of the circuit 80 

Referring now to FIG. 3 the I input to the device is under the control of the microprocessor 70 in accor- 
shown at 50 and the Q input at 51. Each input 50, 51 is 15 dance with software stored in the memory 71. These 
connected to a respective channel sample and hold operations comprise analogue to digital conversion. A 
amplifier 52, 53. The output of the amplifier 52 is con- preselected number of samples can be digitised and 
nected by a signal amplifier and buffer circuit 54 to a processed to obtain statistics data by the microproces- 
display multiplexer and driver circuit 56. The circuit 56 sor 70 in accordance with software stored in the mem- 
is connected to the X deflection circuitry of a cathode 20 ory 71. The statistics can be stored in RAM 72 for sub- 
ray tube 58. The amplifier 53 is connected by a signal sequent further processing. The detail of the calcula- 
amplifier and buffer circuit 60 to a display multiplexer tions which can be carried out by this software will be 
and driver 61. The driver 61 is connected to the Y de- given below. The results of the calculations can be 
flection circuitry of the cathode ray tube 58. Thus it will output for display on the cathode ray tube 58 and can 
be seen that samples from the I channel are displayed 25 also be output to a printer via for example the interface 
along the X axis of the oscilloscope display whilst sam- 82. The interface 82 also permits data from the memory 
pies from the Q channel are displayed along the Y axis. 72 to be output for processing by an external processing 
Each circuit 54, 60 has associated therewith variable device. 

resistors 62, 63 which can be used to adjust the position When the device is set up for a constellation display 
and size of the displayed array. The resistors can be 30 and measurement a graticule can be selected for display 
adjusted by means of knobs on the front panel of the on the cathode ray tube. The operator is provided with 
device shown in FIG. 2. These knobs are shown in the a manual control which allows selection of a graticule 
area 36 on FIG. 2, which is appropriate to the modulation type of the radio 

The output of each amplifier circuit 54, 60 is also transmission. The display of the graticule is controlled 
connected to a re-sampler and comparator 65 which in 35 by the graphics generator 80. In the case of 16 QAM the 
turn is connected by an interface buffer 66 and data bus graticule 32 comprises four equally spaced vertical and 
67 to a microprocessor assembly 68. The microproces- horizontal intersecting lines as shown in FIG. 2. The 
sor assembly has an 8088 device 70 with associated operator is also provided with means for adjusting the 
PROM 71 and RAM 72 with which it can communicate clusters of the constellation relative to the intersections 
via address bus 74 and data bus 75. The data and address 40 of the graticule lines to set the constellation at the ap- 
buses 74, 75 are also coupled to an instrument bus inter- propriate position for a measurement operation. This 
face 78 which is connected to the data bus 67 and also to procedure is carried out manually by the operator using 
an address bus 69. The microprocessor assembly also knobs on the front panel of the device. The knobs adjust 
includes a graphics generator 80 which has outputs to the resistors 62, 63 of circuits 54 and 60 and allow the 
the circuits 56 and 61. An HP-IB interface 82 is also 45 operator to shift the constellation horizontally or verti- 
provided to allow interfacing with other equipment cally or to expand and/or contract the clusters of the 
Other blocks provide clock signals and timing signals constellation to fit the grid arrangement of the graticule, 
for purposes which will be apparent to those skilled in Examples of impairments which may be identified 
the art. from a constellation display are: 

The device also includes a keyboard assembly shown . 50 Demodulator phase lock error (rotation) 
generally at 85. This comprises a keyboard interface Modulator/demodulator quadrature error (orthogonal- 
which is coupled to the data bus 67 and address bus 69 ity error) 

and provides an interface between these buses and man- Interference from an unmodulated (CW) Source, 
ually operable controls on, for example, the front panel Amplitude compression/expansion (e.g. in power am- 
of the device. 55 plifler) 

The timing signals for the sample and hold circuits 52, Amplitude to phase conversion (e.g. in power amplifier) 
53 are derived from clock signals applied to an input 90 Faults or filter misadjustment in one baseband signal 
from where they are fed to a clock signal processing path 

circuit 91 which connects to a time base generator and Examples of constellations which can be produced 
delay circuit 92 and then to multiplexers and latches 94. 60 are shown in FIG. 4. In FIG. 4A the constellation is 
The device also includes display and time base sweep D rotated and this is indicative of a recovered carrier 
to A converters and blanking logic circuitry etc., which phase offset due to off frequency error or to recovery 
is shown schematically at 96. circuit fault. FIG. 4B shows a constellation which has 

In use the signals in the I and Q channels 50, 51 are been compressed. This can result from power amplifier 
sampled by circuits 52 and 53 at instants corresponding 65 fault. FIG. 4C shows an out of square array indicative 
to the instants at which the decision circuitry of a digital of a modulator or demodulator fault. The displays 
radio receiver carries out its sampling operation. The shown in FIG. 4 indicate the types of faults which may 
sampling instants can be varied by means of manually . be present in the radio transmission. They do not how- 
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ever provide accurate quantitative measurements. Ac- upon which the device operates to provide measure- 
curate quantitative measurements can be obtained using ments of various constellation parameters. The values 
the present device as will now be described. are accumulated in internal 16 bit and 32 bit registers. If 
Initially the device is set to its measurement mode by too many samples are accumulated in any group one of 
the operator and the graticule selected which is appro- 5 these registers may overflow. When a new sample 
priate to the modulation scheme. The display 30 of the would cause overflow the device will report an error 
constellation is then adjusted so that the clusters of the discard the sample and terminate the accumulation. The 
constellation correspond closely with the intersections number of samples which can be accumulated before 
of the lines of the graticule. When the device is in- overflow depends upon the input signal and the modu- 
s' 0 ** 1 }° carr y out a measurement operation it ini- io lation scheme. As a worst case, that is when all samples 
tially digitises and processes a large number ofsamples fall m one at ^ extreme measurement limit of 
produced by the samphng cuxuits 52 and 53. The large +312 or _ 512 which ^ extremely unlikely, the limit is 
number can be several hundred or a few thousand. In reached g , 91 sam les> However for a constellation 
g ^l Sui?° with the samples distributed reasonably uniformly 

3 Whence constellation is digitised the samples are 15 7° ng ? * e %T ° U? ?*™t *° ^ * 

, « . ~ . ° . ^Jt larger. In the case of 16 OAM with a reasonable constel- 

selected at random from the contmous stream of input t ~f . . . ^TT " v , * A: T % , 

samples. Two 10 bit analogue to digital conversions ^ lt » ^accumulate on the order of 126000 

detennine the I and Q signal levels of each sample and sam P les ? efore ov J e . r ? 0W oc ? urs " ^ been 

these digital values are stored in the memory 72 of the 20 Previously many digital radio impairments can be dtag- 

microprocessor assembly 68. nosed b ? their effect on constellation display. The 

Referring to FIG. 2 the constellation display area is P resent device ^ lows impairments to be analysed 

square and positioned slightly to the right of the cath- quantitatively and the following is an explanation of 

ode ray tube centre. Signal levels are measured in units ^ee of the quantitative measurements which can be 

scaled to the constellation display area. The centre of 25 made on constellation impairments, 

this display area which corresponds to the centre of the ^ fi«t of these is constellation closure which can 

graticule is the zero point of the I and Q scales. The Q be defined as the RMS size of the clusters divided by 

scale has a range of -512 to +512 from bottom to top one half of the cluster spacing. The computation is 

of the display and similarly the I scale has a range of based upon the statistics described above and stored in 

—512 to +512 from left to right going horizontally. 30 registers in the processor assembly. The computation 

Thus each constellation sample is a pair of numbers can be compared to the N— 1 weighted standard devia- 

which will be in the range —512 to +512. Representing tion of many pocket calculators. In order to obtain a 

the samples as integers facilitates analysis of the constel- single value for the entire constellation all the constella- 

lation. In order to analyse a constellation it is necessary tion clusters are combined together by subtracting the 

to divide the samples up depending on which constella- 35 mean cluster position from the samples in each group, 

tion cluster they belong to. The easiest way to do this is Then the standard deviation is computed using a 

to arrange for the operator to align the constellation weighting of N—M 2 where M 2 is the number of clusters 

relative to the display graticule as described and then ((N— M 2 ) being the number of degrees of freedom). It 

compare the digitised samples to threshold values set will be appreciated that two numbers are obtained for 

midway between each pair of graticule lines. To ^ constellation closure. One is for the I direction and one 

achieve this the measurement scale is based on the posi- for the Q direction. 

tion of each sample on the display rather, than on the The second measurement relates to receiver carrier 

absolute input voltage of each sample. lock angle error which is overall rotation of the constel- 

Thus it will be seen that each digitised sample is alio- lation. The error corresponds to the mean of two phase 

cated to a group whose boundary is defined by four 45 errors m tne digital radio link. These are the orthogo- 

numbers corresponding to the mid points of adjacent nality error between the transmitter I and receiver Q 

graticule lines or the measurement limit in each of the I carriers and the orthogonality error between the trans- 

and Q directions. In the case of 16 QAM there are 16 q ^ reC eiver I carriers. It is usually equal to the 

groups one for each cluster. The digitised samples are error m ^ receiver carrier recovery phase locked 

processed and stored as groups of statistics in the micro- 5Q loop# Anti-clockwise rotation is taken as positive, 

processor assembly for subsequent further processing ^ ^ error is the transmitter/receiver quadrature 

As the samples are digitised the processing assembly Ie error whkh is the ^ ^ ven t0 ^ difference 

carries i out a first processing stage which will be under- from 9Q . of ^ a ]e ^tween the I and Q cluster lines, 

stood by referring to FIG. 5. FIG. 5 illustrates a cluster ^ ds t0 the difference between two phase 

compnsmg a plurality of dots, each dot representing a errors m ^ radio ^ Xhese „ the orth 

sample in the duster For each sample in the cluster ^ error between the transmitter I and receiver Q 

? Vah f * Y r dM r ptCS ! ed X fTV° carriers and the orthogonility error between the trans, 
a co-ordinate system centred on the centre of the dis- ^ Q and receiver jf A ^ mQT 

play. Hie processor assembly calculates for each cius- ^ onds v to ^ mchlded ^ grea ? er ^ 90 - at ^ 
tuL i_ r 1 * 1 * 60 lower left of the constellation. 

The sum of Y 2 constellation statistics previously referred to and stored 

The sum of Xx Y 65 m tDe re S isters ^ described. The computation operates 

These statistics are accumulated in registers, there being ™« mean P?f itio ^ numbe t r «* ^P 1 " 
a separate set of registers for each cluster in the constel- ° f °* c,usteis " ^ «*«nal«w essenUaUy attempts to fit 
lation. These valuw represent the basic statistics data mean cluster P° s,Uons to a model constellation m 
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which the clusters are perfect points, equally spaced in 
both I and Q directions with arbitrary amounts of over- 



8 



OUT and Y OUT is a minimum. In practice the compu- 
tation is a matrix calculation expressed as follows: 



FORI 



A~ x Bj = 



Spacing / 
Cross-modulation / 
Position / 



Ui-d) 2 d - d)Q - d) (/ - <0 

' ill j~\ numbcr of ^P 1 ^^ (/ - d)(/ - d) {f - d) 2 (/ - d) 



M M 



M M 

B/ — X I sum of x values 
i=l;=l 



Ji-d) 



U-d) 



1 



j- 



MODEL 



XOUT 



__YOUT 



25 



30 



all rotation, parallelogram distortion and offset. This 2Q 
model constellation can be described by three vectors 
shown in FIG. 6. The mean squared error values for I 
and Q are defined as the mean over all samples of the 
square of the difference from the ideal constellation 
point The computation determines directly the ideal 
constellation which minimises these error values. FIG. 
7 illustrates schematically this concept. The calculation 
gives the parameters of a model which produces a simu- 
lated output which may be compared with the actual 
output, the parameters being computed so that the least 
mean square error is obtained. From this computation it 
is possible to obtain directly the lock error and quad 
error values. The mathematical expression of these con- 
cepts will be given below. 

If the clusters do not lie on parallel lines so that there 3J 
appear to be four slightly different values for the angle 
between the cluster lines and the graticule lines this is 
not essentially a problem. The ideal constellation al- 
ways has its points on parallel lines so there is only one 
possible value of lock and quad error. The computation ^ 
provides a best fit least mean square ideal constellation 
which gets as close as possible to the measured constel- 
lation. 

Mathematically the model can be considered as fol- 
lows: 



FOR Q 
A~ x Bq = 



Cross-modulation Q 
Spacing Q 
Position Q 



M M 

Bn~1 Z sum of y values (ij) 
/= 1 j= 1 



These evaluations are carried out for all clusters and the 
result is values for the coefficients spacing, cross modu- 
lation and position. The final step is the calculation of 
closure, lock angle and quad angle as follows: 



/CLOSURE . 



/cluster size 
/ spacing 



*2 * 100% 



Lock angle 



45 



Q CLOSURE = P^IT ' 2 • 100% 
Q spacing 

_ . j. _ | (Q cross-modulation) _ 
= L Q spacing 



-1 cross-modulation) 



/ spacing 



50 



Quad angle = —tan' 



.i (O cross-modulation) 
Q spacing 



The values i and j identify the particular transmitter I 
and Q modulation states associated with each constella- 
tion cluster (or sample group). X OUT and Y OUT are 
calculated or predicted values which are compared 55 
with the actual sample values X and Y. X OUT and Y 
OUT can be expressed as follows: 



_ i (/ cross-modulation) 
/ spacing 



OUT = spacing// — d) + cross-modulation X/ — <0" 
-♦-position/ 



60 



OUT = spacing^— ^0+ cross-modulation^/— dy 
+ position^ 

Spacing, cross-modulation and position are coeffici- 
ents which are calculated such that the RMS error 
between predicted X OUT and Y OUT and measured X 



65 



These values when evaluated can be displayed along- 
side the constellation as shown in FIG. 2. The mathe- 
matics explained above is summarised in Appendix A. 

The examples given above relate to the measurement 
of linear effects. It is also possible to use the device to 
measure non-linear effects. Non-linear effects can arise 
from the M-level signals at the modulator, from the 
predistorter and from the travelling wave tube. 

We have identified residual non-linearities which can 
be defined by new measurements. These measurements 
are independent of all linear impairments such as quad- 
rature error and lock error. Three of these measure- 
ments are total residual non-linearity, radial residual 
non-linearity and angular non-linearity. They are 
termed residual because they represent what is left after 
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a least means square error criterion has been used to 
remove all linear effects. They are important because 
the automatic level control circuits of a radio tend to 
compensate for compression and the carrier recovery 
loop tends to compensate for am/pm, any resultant 5 
effect being a residual non-linear displacement of the 
clusters. 

The parameters referred to can be evaluated as fol- 
lows: 

(i) Form means and vacancies for each cluster 10 

(ii) Estimate the state spacing, offset and cross modula- 
tion coefficients 

(iii) Form a modifier matrix from the parameters in (ii) 
above 

(iv) Transform the means in (i) using the modifier ma- 15 
trix to remove the effects of all linear impairments 
The resulting information is shown in FIG. 8. Then 



Q modulation (ij) * 



10 



sum of v values (if) 



- position Q 



number of samples (LJ) 
Q displacement (iJ) = Q modulation (ij) - (spacing Q)(f - d) 



_ ((/modulation a/)) 2 -f (O modulation (i/)) 2 )* 
s(ij) 



' « = («»")- Of ) 

Then the following summations are evaluated over all clusters for 
which (/ — d) or (/ — d) are not both zero. 



m m 
count ~ 2 2 w(ij) 
/=U=1 



Total Residual Non-Lincari 



Radial Residual Non-Linearity 



ity =\f 



2X ■ X 



Number of States 



*RR 



Number of States 



Angular Residual Non-Linearity < 
where 

X « Vector Separation 
R - Radial Error 
A = Angular Error 



2A -A 



Number of States 



35 



40 



It will be appreciated that these parameters can be 
evaluated using the statistics accumulated in the regis- 
ters of the processing assembly. They can be output on 
to the display or to a printer. 

A further set of three non-linear measurements which 
can be made are amplitude compression, amplitude to 
phase conversion, and worst decision margin reduction. 
These all make use of the least mean squared error 
model already described. 

Amplitude compression measures the radial displace- 
ment of the clusters at the four corners of the constella- 
tion as a fraction of the radial distance from the constel- 
lation centre. This is obtained using linear regression of 
the measured radial distance factor u (i j) upon the mod- 
ulation state radial distance s (i j) (see FIG. 9). 

Amplitude to phase conversion measures the angular 
displacement of the clusters at the four corners of the 
constellation. This is obtained using linear regression of 
the measured angular displacement error v (i j) upon the 
modulation state radial distance s (ij) (see FIG. 10). 

Worst decision margin reduction is defined as the 
worst displacement of the mean position of any of the 
clusters from the model position, divided by one half of 
the overall cluster spacing. Two values are computed 
for I and Q and the larger is the worst decision margin 
reduction. 

These measurements are computed using values from 
the least mean square error model, and additional values 
as described below. 

For each cluster, the following are evaluated: 

/«***« «» - r^a - ' 65 

/ displacement (ij) =* / modulation (ij) — (spacing J) (i - d) 



20 sum_s = ^ 2 1 J£^ w(ij) s(ij) 

sum_ss J^Xj u<v) s(ij) s(ij) 

25 mm 

sum u = 22 wO/) u(U) 
i=\j=l J vv ' 

sum_su - Jj^yjj wiij)s(ij) "(ij) 

sum_v - 2 2 w(ij) v(LJ) 
/= 1 1 

sum_sv - .2^ w(ij) s(tj) v(ij) 

There are two alternative choices for the weighting 
factors w(i j). 



30 



50 



60 



number of samples (/./) — 1 



((sura of x 2 ) — (sum of x)V(number of samples)) 

The choice affects the speed of computation and the 
precision of the results. 
Finally the measurements are computed. 

worst decision margin reduction / (WD MR J) = 

maximum / displacement (Li) # 2 , lQ(J% 
spacing / 

worst decision margin reduction Q(WDMR Q) = 



maximum Q displacement ffcft , , 

spacing Q '° 



worst decision margin reduction = 

maximum (WD MR /. WD MR Q) 
sum_u sura ss — sum su sum_s 



am intercept 
am slope = • 
pm intercept 



count sura ss (sum s) 2 

count sum — su — sura u sum— s 

count sum ss (sum— s) 2 

sum_v sum— ss — sum— sv sum— s 
count sum— ss — (sum__s) 2 
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11 

•continued 

co unt sum— sv — sum— v sum— s 
count sum _ss (sum s) 2 



4,918,708 



amplitude compression < 



am slope * (d - 0 * surtjl) 
am intercept 



100% 



amplitude to phase _ p(n ^ , (d _ , } . ^ (2) 
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The above description relates to a device for analy- 
sing 16 QAM transmissions. It will be appreciated that 
it can also be used to analyse many other types of modu- 
lation schemes such as QPSK, 4 QAM, 9 PRS, 49 PRS, 
7 LPR, 64 QAM, 128 QAM, 256 QAM and 512 QAM. 

APPENDIX 



Measurement samples 



Constellation statistics 



(x ft yr)r~ 1 toN 

Samples sliced into buckets 



(xihkf yijid / « 1 to m 
j = 1 to m 
k = 1 to nij 



^Statistics accumulated 
sum-count (y) = n# 



sum-xx(y) = 



sum-yy (*/) = ^2 j j^* 



n-xy 0/) = ^ *jf*7fl* 



^Intermediate calculation 



/ clustor size = 



i T m m [sum-x(y")] 2 "1 



1 ["mm [sum-y(v)] 2 T 

Ocluaersize = vj ^-^^.2, .S 1 sum-yyW) - 8um coumW) J 



m m 

matrix A = ^ ^1 ^ sum-count(y) 



l(t-d) (/-<*) 1 . 



Matrix Bi ss sum-x(y) | (/ - 



matrix 5q = 2^5^ sum-y(y) | (/ - ' 



,[;::] 
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-continued 
Solution of matrix equations 



Matrix 7/ = 



Matrix Tq - 



spacing / 

cross-modulation / 
position / 

cross-modulation Q 
spacing Q 
Postion Q 



^computation of final results 
Constellation analysis remits /eons , e i lalion dosure _ 2 

J spacing 

C constellation closure = 2 ? °'f' er . *» 

y spacing 
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Lock „ g .e = , [«.- ( ^ g- )+ — -' ( / ?SS t )] 



General Notation 
X r I components of sample No r 

Yr Q component of sample no r 35 
n# The number of samples which, after slicing, fell 

within the boundaries of bucket (i j) 
N The total number of samples in the constellation 

measurement. 

m Number of modulator levels in the modulation 40 
scheme, e.g. 4 for 16 QAM. 

d 

»L±± 45 

Xiyjfc I component of sample no k in group (i j) 
Yyjt Q component of sample no k in group (ij) 

NOTES 50 

* Matrix A is positive definite (assuming any three slic- 
ing buckets, not in a straight line, are nonempty) 
therefore the matrix equation always has a solution. 

+ The equations for lock and quad angle are arranged 
to compensate for any difference between I and Q 55 
spacing. The assumption is that the I and Q channels 
have equal modulation level spacing. 

T The use of (i— d) and (j— d) reduces the dynamic 
range, making implementation simpler and more effi- 
cient. 60 

* The equations are solved using L-U decomposition as * 
this improves computational efficiency. 

What is claimed is: 

1. An apparatus for analyzing digital radio transmis- 
sions comprising: 65 

sampling means for sampling received radio signals to 
produce for each sampling instant a signal or sig- 
nals as samples which are representative of the 



modulation state of the transmission at the sam- 
pling instant, and processing means which is con- 
figured to receive and measure a given number of 
said signals, wherein said processing means is con- 
figured to process said signals in accordance with 
one or more stored routines to generate one or 
more parameters which are indicative of the condi- 
tion of the transmission, wherein said samples are 
represented by digitally encoded numbers and said 
processing means is arranged to allocate said num- 
bers to one of a plurality of groups, the number of 
groups corresponding to the number of modulation 
states of the transmission, and wherein said pro- 
cessing means is further configured to carry out a 
plurality of processing stages, a first processing 
stage comprising the evaluation of a set of statistics 
which comprise data representative of the position 
of each sample in a co-ordinate system, and a sec- 
ond stage of which comprises forming, using a least 
mean square error technique, a model cluster array 
which is a closest fit to the array represented by 
said statistics, said modelling stage producing a set 
of coefficients from which is calculated parameters 
representing distortions of a constellation, and 
wherein said processing means is further config- 
ured to evaluate a radial residual non-linearity pa- 
rameter and an angular residual non-linearity pa- 
rameter, said parameters being evaluated by re- 
moving linear effects using a least mean squares 
error technique. 

2. An apparatus for analyzing digital radio transmis- 
sions comprising: 

sampling means for sampling received radio signals to 
produce for each sampling instant a signal or sig- 
nals as samples which are representative of the 
modulation state of the transmission at the sam- 
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pling instant, and processing means which is con- 
figured to receive and measure a given number of 
said signals, wherein said processing means is con- 
figured to process said signals in accordance with 
one or more stored routines to generate one or 5 
more parameters which are indicative of the condi- 
tion of the transmission, wherein said samples are 
represented by digitally encoded numbers and said 
processing means is arranged to allocate said num- 
bers to one of a plurality of groups, the number of 10 
groups corresponding to the number of modulation 
states of the transmission, display means for dis- 
playing said parameters, and wherein said process- 
ing means is further configured to process said 
parameters for display of both diagrammatic and 15 
numerical quantitative information simultaneously 
on said display, wherein said information is indica- 
tive of the condition of the transmission, and 
wherein said processing means is arranged to eval- 
uate the root mean square size of each group of 20 
samples to thereby provide a quantitative measure- 
ment of radio communication impairment, and 
wherein said processing means is arranged to eval- 
uate a radial residual non-linearity parameter and 
an angular residual non-linearity parameter, said 25 
parameters being evaluated by removing linear 
effects using a least mean squares error technique. 

3. An apparatus for analyzing digital radio transmis- 
sions comprising: 

sampling means for sampling received radio signals to 30 
produce for each sampling instant a signal or sig- 
nals as samples which are representative of the 
modulation state of the transmission at the sam- 
pling instant, and processing means which is con- 
figured to receive and measure a given number of 35 
said signals, wherein said processing means is con- 
figured to process said signals in accordance with 
one or more stored routines to generate one or 
more parameters which are indicative of the condi- 
tion of the transmission, wherein , said samples are 40 
represented by digitally encoded numbers and said 
processing means is arranged to allocate said num- 
bers to one of a plurality of groups, the number of 
groups corresponding to the number of modulation 
states of the transmission, display means for dis- 45 
playing said parameters, and wherein said process- 
ing means is further configured to process said 
parameters for display of both diagrammatic and 
numerical quantitative information simultaneously 



60 
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on said display, wherein said information is indica- 
tive of the condition of the transmission, and 
wherein the processing means is arranged to carry 
out a plurality of processing stages, a first process- 
ing stage comprising the evaluation of a set of sta- 
tistics which comprise data representative of the 
position of each sample in a co-ordinate system, 
and a second stage of which comprises forming, 
using a least mean square error technique, a model 
cluster array which is a closest fit to the array 
represented by said statistics, said modelling stage 
producing a set of coefficients from which is calcu- 
lated parameters representing a measurement of 
radio communication impairment, and wherein 
said processing means is arranged to evaluate a 
radial residual non-linearity parameter and an an- 
gular residual non-linearity parameter, said param- 
eters being evaluated by removing linear effects 
using a least mean squares error technique. 

4. An apparatus for analyzing digital radio transmis- 
sions comprising: . 

sampling means for sampling received radio signals to 
produce for each sampling instant a signal or sig- 
nals as samples which are representative of the 
modulation state of the transmission at the sam- 
pling instant, 

processing means which is configured to receive and 
measure a given number of said signals, wherein 
said processing means is configured to process said 
signals in accordance with one or more stored 
routines to generate one or more parameters which 
are indicative of the condition of the transmission, 
wherein said samples are represented by digitally 
encoded numbers and said processing means is 
arranged to allocate said numbers to one of a plu- 
rality of groups, the number of groups correspond- 
ing to the number of modulation states of the trans- 
mission, and wherein said processing means is fur- 
ther configured to evaluate the root mean square 
size of each group of samples to thereby provide a 
quantitative measurement of radio communication 
impairment, and wherein said processing means is 
arranged to evaluate a radial residual non-linearity 
parameter and an angular residual non-linearity 
parameter, said parameters being evaluated by 
removing linear effects using a least mean squares 

error technique. 

***** 
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